Journal of Research of the National Bureau of Standards 

The NBS Scale 
Of Radiance Temperature 



Volume 92 



Number 1 



January-February 1987 



William R, Waters, 

James H* Walker, 

and Albert T. Hattenburg 

National Bureau of Standards 
Gaithersburg, MD 20899" 



This paper describes the measurement 
methods and instrumentation used in 
the realization and transfer of the In- 
ternational Practical Temperature 
Scale (IPTS-68) above the temperature 
of. freezing gold. The determination of 
the ratios of spectral radiance of tung- 
sten-strip lamps to a gold-point black- 
body at a wavelength of 654.6 nm is 
detailed. The response linearity , spec- 
tral responsivity, scattering error, and 
polarization properties of the instru- 
mentation are described. The analysis 



of sources of error and estimates of 
uncertainty are presented. The as- 
signed uncertainties (three standard 
deviations) in radiance temperature 
range from ±2 K at 2573 K to ±0.5 
K at 1073 K. 
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1. Introduction 



Temperatures above the freezing point of gold 
(1337.58 K) are defined on the International Practi- 
cal Temperature Scale (IPTS-68) [l] 1 in terms of 
the ratio of spectral radiances of two blackbody 
sources, one of which is maintained at the gold 
point. Thus, the accurate measurement of the spec- 
tral radiance ratio of two sources is required for the 
realization and dissemination of the radiance tem- 
perature scale. Currently, such measurements are 
performed at NBS using a spectroradiometer oper- 
ated by a computer-controlled data acquisition sys- 
tem, which permits rapid comparison of the 
spectral radiances from a variety of sources. With 
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this system, lamp standards of radiance tempera- 
ture are calibrated in the range 1337,58 K to 2573 
K. For convenience, the calibrations are extended 
to 1073 K, where the scale is defined in terms of 
thermocouples. Calibration services are also avail- 
able for optical pyrometers. The purpose of this 
paper is to provide a detailed description of the 
instrumentation and procedures used to realize and 
disseminate the radiance temperature scale. De- 
tailed information on the calibration services, ap- 
pears in a separate report [2]. 



2. Basic Theory 

Spectral radiance is the radiant power contained 
in a defined area, solid angle, and wavelength inter- 
val, 

L x =d*<P/dA *cose*dtt.dX (1) 

where L x is the spectral radiance, 4> is the radiant 
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flux, A is the area, 6 is the angle between the sur- 
face normal and the direction of propagation, ft is 
the solid angle about that direction, and X is the 
wavelength. The relation between spectral radi- 
ance, wavelength and temperature is given by 
Planck's Law, 



U=c l /ir-\ 5 .[exp(c 2 /\-T)-i\ 



(2) 



where C\ and c 2 are the first and second radiation 
constants, \ is the wavelength in vacuo, and T is 
the temperature. The defining equation for the 
IPTS above 1337,58 K is therefore 

r=L K (T)/U(T Au ) 

= [exp(c 2 AT Au ) - l]/[exp(c 2 A.D - 1] (3) 

where Lx(T) and L\(T Au ) are the spectral radiances 
of the two blackbodies at temperatures T and T Au , 
r Au is the temperature of freezing gold defined as 
1337.58 K, r is their ratio, and c 2 is defined as 
14388 cm-K. In principle, a measurement of the 
ratio at a discrete wavelength with a linear re- 
sponse instrument yields the value of T. 

In practice, the radiance temperature scale is re- 
alized with an instrument of finite bandpass, and an 
integral form of eq (2) must be used. 



r= fL K (TyR y dX/JU(T A ^R k -dX 



(4) 



where R k is the spectral responsivity of the instru- 
ment and includes the spectral transmittance of the 
wavelength-limiting element (e.g., interference fil- 
ter or monochromator), the spectral transmittance 
of all other optical elements, and the spectral re- 
sponsivity of the detector. A determination of R x 
allows a calculation of T from eqs (2 and 4). Deter- 



mination of the ratio for a range of values provides 
a temperature scale over the corresponding tem- 
perature range. 

The radiance temperature scale is typically main- 
tained and disseminated on tungsten-strip lamps, 
which possess a repeatable current vs. radiance- 
temperature relationship not available in present 
variable-temperature blackbodies. Realization of 
the scale with these lamps and a gold-point black- 
body is also practical, and is the typical procedure 
followed by standards laboratories [3]. This scale is 
valid only for the wavelength of realization, since 
the spectral distributions of the lamps are not 
known functions, Traditionally the wavelength of 
realization has been near 650 nm, a convenient re- 
gion for visual optical pyrometers. At NBS, a spec- 
troradiometer system is presently being used and a 
wavelength of 6546 nm (wavelength of a thorium 
spectral line) has been chosen. This method re- 
quires that the instrument relative spectral respon- 
sivity function extend only over an acceptably 
small spectral range, or is known accurately 
enough to determine the wavelength at which the 
integrands of eq (4) have the same ratio as the inte- 
grals. 



3. Measurement Apparatus 

The radiance temperature calibrations described 
in this report are performed on the NBS Facility 
for Automated Spectroradiometric Calibrations 
(FASCAL), This system employs a prism-grating 
double monochromator, mirror optics, a highly sta- 
ble external tungsten reference source, and a photo- 
multiplier tube operated in a linear response mode. 
A block diagram of the system is shown in figure 1, 
Prior calibrations were carried out on the NBS 
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Figure 1-FASCAL radiance temperature 
measurement. 
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Photoelectric Pyrometer [4]. A new version of this 
instrument is under development, which will em- 
ploy interference filters (1 nm bandpass) and re- 
fractive optics. The measurement process will be 
nearly identical to that of the FASCAL system, 
utilizing an external reference source and a linear- 
response photomultiplier tube. A schematic dia- 
gram of the new instrument is shown in figure 2. 

3,1 Gold-Point Blackbody 

The Gold-Point Blackbody (GPBB) is a graphite 
cylinder with a small viewing hole (diameter 1 
mm) at one end, and a conical cavity at the other. 
The cylinder is surrounded by 0.99999 pure gold, 
and the crucible containing the gold is surrounded 
by heating coils and an insulated case. The con- 
struction and characterization of the GPBB have 
been detailed in prior reports [4]. The estimated 
emissivity is 0.9999. The duration of a melt or 
freeze plateau is approximately 5 minutes, and the 
time delay between these observation periods is 
about 5 minutes. 



3.2 Lamp Sources 

Tungsten-strip filament lamps are used both in 
the realization of the radiance temperature scale 
and as secondary standards for scale dissemination. 
Each filament has a small notch in one edge, about 
midway along its length, to aid in determination of 
the filament portion to be calibrated (target area). 
A small mark is placed on the rear of the lamp 



Obj. lens 



Field 
stop 



envelope to permit reproducible angular position- 
ing. 

3.2.1 Vacuum Lamps: Vacuum tungsten-strip 
lamps of the Quinn-Lee type [5,6] are used in the 
scale realization. One lamp is operated at a single 
current to produce a spectral radiance equal to that 
of the gold-point blackbody at 654.6 nm. A second 
lamp is operated at a single current to produce a 
spectral radiance of about eight times that of the 
gold-point blackbody at 654.6 nm (about 1530 K 
radiance temperature). Both lamps are stable to 
better than 0.02% over 200 burning hours when 
operated under these single-level conditions. 

3.2.2 Gas-Filled Lamps: Gas-filled tungsten- 
strip lamps are calibrated as secondary standards 
for dissemination of the radiance temperature scale. 
Prior to calibration, the lamps are annealed at a 
radiance temperature of 2350 °C for 2 hours on di- 
rect current. Each lamp is examined for the varia- 
tions in spectral radiance with angle of emission, 
while set at a radiance temperature of 1700 °C. The 
lamp orientation is chosen to minimize variations in 
spectral radiance with rotation, and an alignment 
arrow is then etched on the rear of the envelope. 
The amount and orientation of polarization effects 
are also determined. Lamps selected for calibration 
display a degree of polarization of 0.005 or less. 

3.3 Spectroradiometer 

3.3.1 Fore-optics: Sources are imaged with unit 
magnification at the entrance slit of the monochro- 
mator by the two front-surface aluminized mirrors 
shown in figure 1. The mirror surfaces are stripped 
and recoated at intervals to reduce signal loss due 
to oxidation. The plane mirror directs the beam to 
the spherical mirror (radius of curvature 1220 mm) 
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Figure 2-Schematic of photoelec- 
tric pyrometer. 
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along a line about 1.5 degrees from the spherical 
mirror axis. The spherical mirror focuses the beam 
onto a mirror-surfaced mask, placed immediately in 
front of the entrance slit. The mask is engraved 
with horizontal and vertical scales with 0. 1 mm di- 
visions, and is viewed at high magnification by ei- 
ther a telescope or a video camera to allow for 
precise positioning of the sources. The mask deter- 
mines the height of the system field stop (source 
target area). The entrance slit determines the 
width. The stop dimensions are 0,8 mm high by 0,6 
mm wide for the scale realization and transfer. The 
system aperture stop is located within the 
monochromator. 

3.3.2 Monochromator: A prism-grating double 
monochromator is employed to minimize spectral 
scattering and avoid multiple orders. The disper- 
sion is about 4 nm/mm at the 654.6 nm wavelength 
setting. The entrance aperture (solid angle) is 
rectangular in shape, with a vertical angle of 0.125 
radians and a horizontal angle of 0.0625 radians. 
The wavelength setting is calibrated against a spec- 
tral line standard (thorium) $nd is repeatable to 
within 0.05 nm. The entrance, intermediate, and 
exit slits are adjustable from 0.01 to 3.0 mm as a 
unit, resulting in a nearly triangular-shaped spectral 
responsivity function. 

3*3,3 Detector: For the radiance temperature 
determinations an end-on 11 -stage photomultiplier 
with quartz window and S-20 spectral response is 
placed behind the exit slit. The detector is cooled 
to 258 K with a thermoelectric cooler. The anode 
current is amplified and converted to a 0^10 volt 
signal by a programmable amplifier. 

3.4 Control and Data Acquisition System 

The FASCAL system employed for the radiance 
temperature calibrations permits control of the en- 
tire measurement process from a remote operator 
console after initial source alignment. Component 
positions, instrument settings, sequence of opera- 
tions, and data collection are effected by either 
stored computer programs, operator commands, or 
a combination of the two. 

The system is directed by a microcomputer and a 
high-speed disk system for program and data stor- 
age. A modular interface controller [7] provides 
the link between instruments and computer. All 
measurement signals are multiplexed into the digi- 
tal voltmeter through the interface scanner, and the 
instruments are remotely programmed and con- 
trolled through interface modules. All instrument 
settings and signal outputs are printed and stored 
on disk for later analysis. The spectroradiometer 



(fore-optics, monochromator, and detector) and a 
closed-circuit TV camera are mounted on a car- 
riage. The carriage can be positioned by remote 
command along a linear track, to align the spectro- 
radiometer with one of the sources mounted at 
fixed stations along the track. The average move 
time between stations is a few seconds, and posi- 
tions are repeatable to about 0. 1 mm. The TV cam- 
era presents a highly magnified image of the 
monochromator entrance slit mask to video dis- 
plays at the spectroradiometer and at the operator 
console for initial source alignment and subsequent 
monitoring. 

4. Measurement of Instrument and Source 

Parameters 
4.1 Spectral Responsivity Function 

The relative spectral responsivity function of the 
spectroradiometer is determined by the indirect 
method [8]. In this method, the relative responsiv- 
ity function is treated as the product of two terms, 
the responsivity factor and the slit-scattering func- 
tion, where the responsivity factor depends only 
upon the wavelength of the observed flux and the 
slit-scattering function depends only upon the dif- 
ference between the wavelength setting of the 
monochromator and the wavelength of the flux. 
This factorization of the spectral responsivity func- 
tion is valid if the instrument dispersion, aberra- 
tions, scattering, and diffraction are constant over 
the wavelength region of interest. This assumption 
is valid in the central portion of the relative respon- 
sivity function, but values for the distant wings are 
subject to error due primarily to changes in scatter- 
ing and dispersion. 

The responsivity factor is obtained by spectrally 
scanning a continuous source standard of spectral 
radiance with narrow (0.1 mm) slits. To determine 
the slit-scattering function, an integrating sphere ir- 
radiated by a krypton or argon laser is spectrally 
scanned by the spectroradiometer, with the slit 
widths set at the 0.6 mm width used in the scale 
realization and transfer. The plot of the output sig- 
nal versus wavelength is the mirror image of the 
plot of the slit-scattering function versus wave- 
length [8], For a 647 nm krypton laser, the function 
is nearly triangular in shape with a width at half- 
height of 2.2 nm. Relative to the peak value, the 
measured values decrease to about 10™ 3 at 3 nm, 
10 4 at 15 nm, and 10 7 at 70 nm from the central 
wavelength. At 150 nm from the central wave- 
length, the value decreases to 10™ 8 in the short- 
wavelength wing and to 10~~ 9 in the long- 
wavelength wing. Scans with 488 and 514 nm (ar- 
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gon), and 676 nm (krypton) yield similar results. 
These values were confirmed over the central and 
near wings portion of the function by measure- 
ments with the direct method, using a dye laser 
tuned through a series of wavelengths with the 
spectroradiometer set at a fixed wavelength [9]. 
Since the function changes very slowly with wave- 
length in the visible region, the measurement at 647 
nm yields the slit -scattering function at 654.6 nm. 

4.2 Linearity of Response 

The degree of linearity of the spectroradiometer 
response is determined with an automated beam 
conjoiner [10]. A beam from a constant source is 
split into two branches, whose fluxes are indepen- 
dently attenuated or blocked before recombination 
and further attenuation. The flux contribution from 
both branches is equal to the sum of the fluxes from 
each branch when measured separately (additiv- 
ity). The device provides 96 levels of flux over a 
range of about 1 to 500. The levels are presented in 
random order to avoid systematic errors, and are 
interspersed with 29 zero flux levels. A microcom- 
puter controls the attenuating filters and records 
the filter positions and radiometer signals. The data 
is least-squares fitted to a polynomial response 
function, to determine a correction factor by which 
the radiometer output signal must be multiplied to 
obtain a quantity proportional to radiant flux. 

The measured instrument response is linear to 
within ±0.2% for a range of photomultiplier an- 
ode currents from 1 to 500 nanoamperes. For cur- 
rents much less than 1 nanoampere, the signal is 
limited by noise. For currents greater than 1 mi- 
croampere, the linearity correction increases 
rapidly, rising to 3% at 7 microamperes. The anode 
current is restricted to less than 500 nanoamperes 
during measurements by selection of appropriate 
photomultiplier tube voltage. Correction factors 
for the amplifier ranges are determined from mea- 
surement of a known electrical current and com- 
bined with the linearity correction factor. 



measure the polarization properties of the lamp 
sources. 

The determination of spectroradiometer and 
sheet polarizer properties consists of spectral radi- 
ance measurements of the sphere source alone, 
with a pair of similar polarizers interposed in the 
beam and set at a number of angular positions (ro- 
tations about the optic axis), and with each polar- 
izer individually placed in the beam at the same 
angular positions. In order to account for departure 
from ideal behavior, a Mueller transmittance ma- 
trix [11] containing six parameters is assumed for 
the sheet polarizers. Circular polarization is as- 
sumed to be negligible. The spectroradiometer po- 
larization direction is determined in a preliminary 
experiment and chosen as the polarization refer- 
ence direction, leaving only the degree of polariza- 
tion to be determined for the instrument. Our 
measurements provide about 200 equations involv- 
ing 10 unknowns, whose values are then obtained 
by a non-linear least-squares solution. 

The sphere source is replaced by the lamp source 
whose properties are to be measured using a char- 
acterized sheet polarizer and spectroradiometer. 
Measurements are made with the lamp alone and 
with the polarizer set at the same angular positions 
as before. This results in 25 equations involving the 
two source unknowns, whose values are obtained 
by a least-squares solution. The polarization of the 
lamp source is specified by the direction of maxi- 
mum linear polarization r, and the degree of polar- 
ization 



■* — f^max ^mirJ/vA'max -\-J^ m i n ) 



(5) 



where £ max and L min are the maximum and mini- 
mum readings of a polarization-indifferent ra- 
diometer when an ideal linear polarizer is rotated 
in the source beam. The factor (see Appendix) re- 
quired to reduce the signal ratio of eq (4) to a value 
which would be obtained with a polarization-indif- 
ferent radiometer is 

(1 +£/>.)-' 



4.3 Polarization 

The polarization properties of the spectrora- 
diometer and the gas-filled lamps are measured 
with the aid of dichroic (linear) polarizers posi- 
tioned in motorized rotating mounts. The sheet 
polarizer properties and those of the spectrora- 
diometer are determined in an initial set of experi- 
ments, using an illuminated integrating sphere as a 
source of unpolarized radiation. The characterized 
polarizer and spectroradiometer are then used to 



where />i = p-cos2r and B is the spectroradiometer 
polarizance (degree of linear polarization intro- 
duced by the spectroradiometer). For the spectro- 
radiometer employed here, the measured 
polarizance is 0.26 at 654.6 nm. The degree of po- 
larization of a typical lamp selected for calibration 
is about 0.003. The uncertainty in the radiometer 
polarizance B is estimated as ±0.0005, and the 
lamp polarization uncertainties are estimated as 
±0.002 (uncertainties are stated at the three stan- 
dard deviations level). 
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4.4 Size of Source 

The "size of source" effect (signal contribution 
due to flux which originates outside the target area 
and is scattered into the measured beam by the 
fore-optics) is determined by observing the signal 
from a 0.6 by 0.8 mm segment of a uniform diffuse 
source, and noting the change in signal when the 
surrounding area of the source is changed by plac- 
ing various masks on the diffuse source. The masks 
expose source areas which closely approximate the 
radiant areas of the lamp and blackbody sources 
used in the scale realization and transfer. As a 
check, the effect is also evaluated by observing 
changes in the near-zero signal from a "black hole" 
(an absorbing cavity slightly larger than the 0,6 by 
0.8 mm field stop) as the various surrounding area 
masks are positioned. The observed differences are 
used to apply a correction to the signals observed 
in source comparisons. The measured effect varies 
from 0.04% to 0.1% at 654.6 nm depending upon 
the elapsed time since thejast mirror r.ecoating. 



5. Process of Realization and Transfer 

The process of realization and transfer consists 
of three steps, all carried out at a wavelength of 
654.6 nm. In step one, the gold-point spectral radi- 
ance is transferred to a Quinn-Lee vacuum lamp. In 
step two, a second vacuum lamp at about 1530 K 
radiance temperature (about eight times the gold- 
point spectral radiance) is compared to the first 
vacuum lamp. In step three, the second lamp is 
compared to the test lamp. After resetting the test 
lamp current, step three is repeated for each de- 
sired radiance temperature. If the vacuum lamps 
remain constant between observations, the product 
of the signal ratios from the three comparisons is 
the signal ratio of the test lamp to the GPBB. With 
appropriate corrections for size of source and po- 
larization, the spectral radiance and thus the radi- 
ance temperature of the test lamp can be evaluated 
from eqs (4 and 2). The polarization and the spec- 
tral radiance distributions over ft, A , and X of the 
vacuum lamps are of no concern here, since the 
lamps serve only to reproduce the spectroradiome- 
ter signal between comparisons. The only require- 
ment is that the vacuum lamp parameters remain 
constant between comparisons with the gold-point 
blackbody (GPBB) and with the test lamp. This is 
satisfied by the constant-current operation of the 
stable (0.02% per 200 hours) vacuum lamps. Com- 
parison with the GPBB is performed about every 



50 to 100 hours of lamp operation, and the vacuum 
lamps are compared with each other weekly dur- 
ing experiments. The use of the first vacuum lamp 
avoids the inconvenience of manipulating the 
GPBB through its melting and freezing point cy- 
cles and provides a continuous measure of gold- 
point radiance. It also increases the number of 
experiments allowed between maintenance or re- 
placement of the gold-point furnace. The use of the 
second vacuum lamp keeps all signal ratios within 
the linear response region of the spectroradiometer 
(i.e., within the t to 500 nanoampere range of an- 
ode current). 



6, Data Analysis and Uncertainties 

6.1 Temperature Values 

The signal ratio of the test lamp to the GPBB is 
obtained from the product of the signal ratios mea- 
sured in the first vacuum lamp comparison with the 
GPBB, the second (1530 K) vacuum lamp com- 
parison to the first, and the test lamp comparisons 
to the second lamp. The ratio is then multiplied by 
correction factors to account for the size-of-source 
disparity, polarization error, and departure from 
linear response where appropriate. An estimated 
test lamp radiance temperature is then calculated 
from the ratio of Planck functions for two black- 
bodies which have the same signal ratio (one black- 
body at the gold point), using 654.6 nm as the 
wavelength. This temperature can then be used to 
determine by iteration the exact temperature which 
will satisfy eq (4). For the spectral responsivity 
function of this spectroradiometer, this exact tem- 
perature differs from the estimated temperature ob- 
tained from the Planck ratio by an amount which is 
small (less than 0,2 K) and is a simple function of 
the temperature. Therefore, to avoid the repetitive 
iteration process, the temperature calculated from 
the Planck ratio is corrected to the desired eq (4) 
value by this known difference. 

6.2 Uncertainties 

The uncertainties in the radiance temperature 
values assigned to the calibrated lamps are obtained 
from the observed imprecision of the measure- 
ments and the estimated systematic error in both 
the measured and the provided quantities (e.g., 
temperature of melting gold). Uncertainties ob- 
tained from observed imprecision and from pub- 
lished values for the physical constants are based 
upon three standard deviations. Uncertainties of 
systematic errors are estimated at the equivalent of 
three standard deviations. 
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In order to examine the contributions of the var- 
ious errors to the uncertainty in radiance tempera- 
ture, an approximate equation for the complete 
measurement process can be developed by using 
the Wien approximation to eq (2), 



L x ^(ci/7r)-X " 5 *exp( - Cz/\' T) 



(6) 



to express the spectral radiance of a blackbody. 
With this approximation, eq (3) becomes 



r = exp(c 2 A- r Au )/exp(c 2 A- T) 
=exp[(c 2 A).(l/r Au ^l/r)]. 



(7) 



Solving for T, and expressing r in terms of the mea- 
sured ratios and their correction factors, we can 
express the complete measurement process as 

r={(i/r Au ) 

-(\/c 2 )Ms-d<f>M -MyM 2 /(l +B.pi)]}' 1 (8) 

where the definitions of the quantities and their es- 
timated 3 or uncertainties are: 
Mo, signal ratio GPBB vs. first vacuum lamp 
(0.12%) 

M h signal ratio first vacuum lamp vs. 1530 K lamp 
(0.12%) 

M h signal ratio 1530 K lamp to test lamp (0.2™ 
0.5%) 

s, size-of source correction for GPBB vs. test lamp 
(0,1%) 

d, correction for test lamp drift during calibration 
(0.1%) 

/, linearity-range factor correction (0.04-0.1%) 
T Au> IPTS-68 temperature of freezing gold (0.4 K) 
c h second radiation constant (0.00014 cm-K) 
J5, spectroradiometer polarizance (0.05%) 
p h lamp polarization component (0,2%) 
X, wavelength setting at 654,6 nm (0.15 nm) 

Radiance temperature uncertainties due to the 
factors of eq (8) are obtained from the partial 
derivative with respect to those factors and the es- 
timated uncertainty in the factor (propagation of 
error). For example, the calculated uncertainty in 
radiance temperature at the 2300 °C (2573.15 K) 
point due to the 0.4 K uncertainty in T Au is 

Ar=(2573.15/1337.58) 2 <0.4)=1.48K. 

Differences between errors calculated by eq (8) 
and those calculated by the exact Planck equation 
are negligible. Note that for the wavelength X this 
process yields the error due to inserting the wrong 



wavelength in the calculation, not the error due to 
incorrect wavelength setting. 

In addition to the factors which appear explicitly 
in this relation, uncertainties in the ratios M Qt M u 
and M 2 arise from errors in the wavelength setting 
X, in the current measurements of the vacuum (0.2 
ma) and gas-filled (2 ma) lamps, in the alignment 
of lamps and in the measured spectral responsivity 
function. The uncertainties in the ratios due to 
wavelength setting and current are assessed by 
measurement of the change in signal ratio when 
varying these quantities. The effect upon the signal 
ratios due to the uncertainty in the measured spec- 
tral responsivity function is determined by solving 
eq (4) for a range of R (A > X) values, using the 
known spectral radiance distribution of the GPBB 
and an approximate test lamp distribution derived 
from spectral scans of such lamps. The radiance 
temperature uncertainties due to these factors are 
then deduced from the ratio uncertainties as before. 
The uncertainties in signal ratio, wavelength set- 
ting, lamp currents and lamp alignment are consid- 
ered random errors; the remaining errors are 
systematics. Table 1 summarizes the uncertainties 
obtained by this process. 
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Table 1, Summary of estimated uncertainties in degrees C. 



Source 

of 
Uncertainty 






Temperature (*C) 






800 


1100 


1400 


1800 


2300 


Signal Ratio M (r) 


0.06 


0.10 


0.15 


0.23 


0.36 


Signal Ratio M { (r) 


0.06 


0.10 


0.15 


0.23 


0.36 


Signal Ratio Mi (r) 


0.27 


0.27 


0.28 


0.37 


0.67 


Size of Source (s) 


0.05 


0.09 


0.13 


0.20 


0.30 


Lamp Drift (s) 


0.05 


0.09 


0.13 


0.20 


0,30 


Linearity (s) 


0.05 


0.03 


0.05 


0.13 


0.30 


Temperature Freezing Gold (s) 


0.26 


0.42 


0.62 


0.96 


1.48 


Second Radiation Constant (s) 


0.02 


0.0 


0.04 


0.11 


0.22 


Polarization (s) 


0.04 


0.06 


0.09 


0.14 


0.21 


Wavelength Setting (r) 


0.08 


0.05 


0.08 


0.20 


0.44 


Vacuum Lamp Current (r) 


0.05 


0.08 


0.13 


0.18 


0.29 


Test Lamp Current (r) 


0.29 


0.18 


0.12 


0.09 


0.07 


Spectral Responsivity (s) 


0.03 


0.01 


0.0 


0.02 


0.03 


Lamp Alignment (r) 


0.13 


0.21 


0.32 


0.49 


0.76 


Total estimated 3<r uncertainty 


0.5 


0.6 


0.8 


1.3 


2.0 


on Thermodynamic Scale 












(square root of sum of squares) 













Note: Random errors denoted by (r), systematic by (s). 



APPENDIX 



Measurement of Instrument and Source Polarization 



The measurement of source polarization and the 
determination of its effect on the measurement of 
spectral radiance is carried out with a sheet polar- 
izer mounted in a motorized rotating mount. This 
mount fits a slot milled into a bracket attached to 
the front of the spectroradiometer. By this means 
the polarizer can be reproducibly positioned in the 
same location between the source and the spectro- 
radiometer whenever required. 

Before the polarizer can be used for source po- 
larization measurements its polarization and that of 
the spectroradiometer must be measured. This is 
accomplished with the help of an integrating 



sphere and a second polarizer similar to the first. 
The integrating sphere is mounted in the radiome- 
ter source plane so that its exit port serves as a 
source of unpolarized radiance. Its entrance port, 
located at about 90° from the exit port, is irradiated 
with a suitable light source, such as a quartz-halo* 
gen lamp. Then a series of spectral radiance mea- 
surements of this source is recorded: a) with no 
polarizer, b) with the first polarizer in its slot, 
c) with the second polarizer added in a second 
milled slot in the mounting bracket between the 
first polarizer and the source, and, d) with the first 
polarizer removed leaving only the second polar- 
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izer. At each wavelength of interest measurements 
are taken as a function of polarizer orientation as 
the polarizers are rotated about the optic axis. The 
first polarizer is typically sampled every 15° and 
the second every 30°, The second polarizer is tilted 
in its mount by about 10° to avoid polarizer-polar- 
izer interreflections. All of these measurements are 
automated and require typically about an hour at 
each wavelength. 

In the analysis of this data we write the first row 
Mueller matrix elements of the spectroradiometer 
responsivity as 



R^Ro 



IB00 



(la) 



where R is the ordinary responsivity neglecting 
polarization and we take the reference direction for 
polarization to be in whatever direction the ra- 
diometer polarizance lies, thus forcing the 45° com- 
ponent to vanish. The reference direction was 
determined to be horizontal by simply observing 
the effect of a hand-held polarizer of known polar- 
ization axis upon a lamp measurement. The polar- 
ization axis of the polarizer, in turn, was confirmed 
by viewing sunlight reflected at glancing incidence 
through the polarizer. Circular polarization is as- 
sumed to be negligible for the type of instrument 
and sources employed here. The quantity B repre- 
sents the degree of linear polarization of the spec- 
troradiometer. The Mueller matrix of a polarizer is 
taken to be of the form 



P = 



s 


dC 


dS 





dQ 


eC 2 +f 


gSC 





dS 


X 


X 


X 





X 


X 


X 



(2a) 



where C=cos 2(<f>™<f> ), S=sin 2(4> — <|> ) and the 
parameters s f d, e t f, g> and <j> are obtained by 
least-squares fitting of the measurements. The an- 
gle <f> is the orientation angle of the polarizer read 
from a scale on the polarizer mount. The matrix 
elements marked X are not involved in these mea- 
surements. This form of the Mueller matrix was 
found, by trial and error, to fit our data well The 
simpler model of an ideal dichroic polarizer in 
which e=s and g=s—f (where, in terms of the 
principle transmittances ki and k 2 , we have 
s = pc 1 +k 2 ]/2 and d = [k ] — k 2 ]/2) is adequate for 
our visible and infrared polarizers but fails to fit the 
ultraviolet polarizer data. Finally, the Stokes vec- 
tor of the unpolarized source is given by 



Ju —~JUq 



The various measurements then lead to equations 
of the following forms. 

No polarizer: V=R >L =RoL (4a) 

Polarizer #1: V($)=R-P>L 

=R L ( fs+BdC) (5a) 

Both polarizers: V($,ty)=R>P>P'>L 

=R L ($s'+Bs'dC+dd f CC+Bd'fC 

+dd'S$'+Bd'gSCS'+Bd'eC 2 C) (6a) 

Polarizer #2: Vty^R-P'-L 

=RoL (s'+Rd'C) (7a) 

where C'=cos 2(<f>'-<V). S' = sin 2(<f>'-(f)o'), and V 
represents the recorded radiometer signal output as 
a function, in general, of polarizer orientation an- 
gles. Altogether our measurements produce about 
200 equations. The values of the ten unknowns s, 
s\ d, d\ e,f t g, B, 4> 0j and <f> ' are then obtained by 
a non-linear least squares solution. As far as polar- 
ization is concerned the parameters s, d, e> /, g, 
and <f> completely characterize the first polarizer, 
at least within the accuracy of the model repre- 
sented by the Mueller matrix of eq (10), and the 
parameter B characterizes the spectroradiometer. 
Since the factor RqL is common to all equations 
these measurements' yield tio information about it. 
Using the characterized polarizer and radiometer 
we are now able to measure the polarization of an 
unknown source by recording spectroradiometer 
measurements obtained from the source: 
a) without a polarizer, and b) with the first polar- 
izer again mounted on the bracket as when it was 
characterized. As before, the measurements are 
made at 15° intervals of the polarizer. We take the 
Stokes vector of the unknown source to be 



La — JLo 



(8a) 



where, as always, we have neglected circular po- 
larization. p x and pi are given by p y =p -cos2t and 
p 2 =p.$in2r> where r is the direction of polarization 
measured from the reference direction of the ra- 
diometer and p is the degree of polarization, 



(3a) 



p — (L max ™L min)/(L max +L min). 



(9a) 
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In eq (17) Lmax and Lmin are the maximum and 
minimum readings of a polarization-indifferent ra- 
diometer when an ideal linear polarizer is rotated 
in the source beam. Then our measurements lead to 
a no-polarizer equation: 



and to 24 equations: 



(10a) 



V(T)=RM(s+BdC) + (dC+BeC l +Bf)p x 
+ (dS+BgSC)p 2 ). 



(Ha) 



Aside from the factor RqL which is common to all 
equations and cancels out, the only unknowns in 
these 25 equations are p x and p 2 which we can then 
obtain by a least-squares solution. These two 
parameters completely characterize the state of 
(linear) polarization of the source. Normally, only 
the second described experiment of 25 measure- 
ments need be performed to measure the polariza- 
tion of a source since the polarizers and radiometer 
are relatively stable and, once characterized, will 
only infrequently need to be measured. 

Equation (10a) gives the relationship between 
the spectral radiance of a source, L , and the ra- 
diometer output signal, V, assuming linearity of re- 
sponse and only linear polarization. If we apply 
this equation both to the measurement of an un- 



known source and to the measurement of a spectral 
radiance standard used to calibrate the radiometer 
we obtain: 



L =mV/VW + WO +Bpd] 



(12a) 



where Li is the known spectral radiance of the 
standard, V/V s is the measured signal ratio, B is the 
radiometer polarizance, and p{ and p x are the rela- 
tive Stokes components measured as described 
above for the standard and for the unknown 
source. (For a blackbody standard, presumably, 
/>f=0.) This shows that the presence of polarization 
introduces the extra factor 

(1 +2^0/(1+^,) 

into the simple measurement relationship which is 
valid in the absence of source polarization or if the 
radiometer is indifferent to polarization. 

The polarization measurement uncertainty by 
this process can be judged by the statistics of the 
least-squares fitting and by the measurement repro- 
ducibility when different polarizers are used. Both 
indications suggest typical uncertainties (one stan- 
dard deviation) in the visible where the polarizers 
are very good and the noise level very low of 
±0.0003 in the radiometer polarizance B and un- 
certainties of ±0,001 in the source relative Stokes 
components p x and /? 2 . 
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